Dynamical moment of inertia and quadrupole vibrations in rotating nuclei 



R.G. Nazmitdinov 1 ' 2 D. Almehed 3 , F. Donau 4 
1 Max-Planck-Institut fur Physik Komplexer Systeme, D-01187 Dresden, Germany 
2 Bogoliubov Laboratory of Theoretical Physics Joint fnstitute for Nuclear Research, 141980 Dubna, Russia 
3 Department of Physics, UMIST, P.O. Box 88, Manchester M60 1QD, United Kingdom 
4 Institut fur Kern- und Hadronenphysik, FZ Rossendorf, 01314 Dresden, Germany 



o 
o 

(N 
-i— > 
o 

O 

o 



> 

OS 

o 
o 



o 



o 

=5 



X 



The contribution of quantum shape fluctuations to inertial 
properties of rotating nuclei has been analysed within the self- 
consistent one-dimensional cranking oscillator model. It is shown 
that in even-even nuclei the dynamical moment of inertia cal- 
culated in the mean field approximation is equivalent to the 
Thouless-Valatin moment of inertia calculated in the random 
phase approximation if and only if the self-consistent conditions 
for the mean field are fulfilled. 

PACS numbers: 21.60.Ev, 21.60.Jz, 21.10.Re 

A description of rotational states is one of the oldest, yet 
not fully solved, problem in nuclear structure physics. 
While various microscopic models based on the cranking 
approach [|IJ|^] describe reasonably well the kinematical 
moment of inertia — — (dE/d£l)/Q for a finite an- 
gular frequency fi, there is still a systematic deviation of 
the dynamical moment of inertia = —d 2 E/d£l 2 (E 
is the energy in the rotating frame) from the experimen- 
tal data at high spins Since the moments of inertia 
are the benchmarks for microscopic models of collective 
motion in nuclei, the understanding of the source of the 
discrepancy becomes a challenge for a many body theory 
of finite Fermi systems. 

The pairing correlations introduced in nuclear physics 
by Bogoliubov jij] and Bohr, Mottelson and Pines || 
improved the description of the inertial nuclear proper- 
ties, especially in the low spin region It was con- 
jectured that at high spins the Coriolis and centrifugal 
forces break Coopers pairs from the pairing condensate 
and cause the transition from the superfluid to the nor- 
mal (unpaired) fluid @. As a consequence, the rigid body 
(kinematical) moment of inertia should be reached at the 
normal phase. However, even for superdeformed nuclei 
where the pairing condensate is expected to be strongly 
quenched, the moment of inertia is usually lower than 
the rigid body value. 

In fact, the conjecture was based on the analogy between 
moment of inertia of a rotating nucleus and magnetic 
susceptibility of a macroscopic superconductor under a 
magnetic field. Therefore, the conjecture may lose its va- 
lidity for finite Fermi systems like nuclei, while remaining 
correct for an infinite number of particles. For example, 
the deviation from the rigid body value could be par- 
tially explained due to shell effects Q. In present pa- 
per we demonstrate that correlations caused by shape 
oscillations of a system are another important ingredi- 
ent which is missing in all state of the art calculations 



of the moments of inertia for rotating nuclei. We fo- 
cus our analysis upon the dynamical moment of inertia, 
since the jW 

contains more information about differ- 
ent properties of the system due to the obvious relation 

j(2) =j(i) + n d jw/dn. 

The first attempt to take into account the contribution of 
the quadrupole oscillations to the correlation energy and 
to the moment of inertia within the cranking+random 
phase approximation (RPA) approach || suffered from 
an inconsistency between the mean field and the residual 
interaction. In addition, the calculations were done in 
a restricted configuration space (only 3 shells have been 
considered) . For a time the solution of the problem was 
postponed, since there was no practical recipe to calcu- 
late the contribution of the correlation energy, even in a 
restricted configuration space. Using the integral repre- 
sentation method developed recently in |lC[ , a full energy 
Erpa can be calculated in the RPA order KIO] with a 
high accuracy and with minimal numerical efforts. Con- 
sequently, this energy can be used to calculate the dy- 
namical moments of inertia . On the other hand, 
in literature it is stated that the dynamical moment of 
inertia should be equivalent to the Thouless-Valatin mo- 
ment of inertia p2fl . Since the realistic application of the 
Thouless-Valatin theory requires the self-consistent solu- 
tion of the mean field and the RPA equations, until now 
this point is not clarified. 

To understand the role of shape oscillations upon the 
value of the moment of inertia and to calculate the 
Thouless-Valatin moment of inertia in a simple but still 
realistic model, we use a self-consistent cranked triaxial 
harmonic oscillator model as a rotating mean field. The 
self-consistent residual interaction constructed according 
to the recipe |b|] is added to describe the collective exci- 
tations in a rotating system. Since all shells arc mixed, 
we go beyond the approximation used in j^] (for a crank- 
ing harmonic oscillator see also |l4|]). Notice that the 
contribution of the pairing vibrations to the correlation 
energy aside of the one from the shape vibrations is also 
important (see Jl5|,|l6| and references there). However, 
there are some open problems with the choice of the self- 
consistent pairing interaction. Therefore, the combined 
effect of the both types of vibrations is beyond the scope 
of the present investigation and we leave this problem for 
the future. 

The mean field part of the many-body Hamiltonian 
(Routhian) in the rotating frame is given by 
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H = J2(ho 



tin — r^Z/ T 



(i) 



where the single-particle triaxial harmonic oscillator 
Hamiltonian h is aligned along its principal axes and 
reads 



h 



-L/ + ^(^ 2 +^ 2 +^z 2 ). (2) 



The eigenmodes and the total energy of the mean field 
Hamiltonian Eq.(pT) are well known p7Hl9|], 



emf = h fo^y^ + ^+x + + ^-X] 



(3) 



(4) 



Here, J2k = (J2j ( n k + V 2 )j) and n k = o-t a k ( k 
:r . + 7 — ) where a£, a% are the oscillator quanta operators. 
The lowest levels are filled from the bottom, which give 
the ground state energy in the rotating frame. The Pauli 
principle is taken into account such that two particles 
occupy one level. The minimisation of the total energy 
Eq.(Q) with respect to all three frequencies, subject to 



the volume conservation condition 



U! x LOyLO z 



yields 



the self-consistent condition |20|,|2l|] for a finite rotational 
frequency 



(5) 



It should be pointed out that the condition Eq.(||) pro- 
vides generally the absolute minima in comparison with 
the local minima obtained from the condition of the 
isotropic velocity distribution ]lq,|l9|| 



(6) 



at large rotational frequency. 

To analyse the contribution of the quadrupole shape os- 
cillations we add to the mean field Hamiltonian Eq.(jl]) 
the self-consistent interaction resulting from small angu- 
lar rotations around the x—, y—, z— axes and small vari- 
ations of two the intrinsic shape parameters e and 7 [ p"3| . 
Consequently, the total Hamiltonian can be expressed as 



H-rpa — Hq — HL X — — 



2 

E 



H - fiL 



(7) 



The effective interaction restores the rotational invari- 
ance of the Hamiltonian Hq such that now [H, Li] = 
(i = x, y, z) in the RPA order. The self-consistency 
condition Eq.(^|) fixes the quadrupole strength k — 



T?--7pry-. Here a mean value 



(x + y ' + z' 



and quadrupole operators — r 1 Y2 ti are expressed 



in terms of the double-stretched coordinates qi = ^-q%, 
(qi — x, y, z). We remind that the self-consistent residual 
interaction does not affect the equilibrium deformation 
obtained from the minimisation procedure. 
Using the transformation from pi , variables to the 
quanta , a k Jl9[ ] , all matrix elements are calculated 
analytically. We solve the RPA equation of motion for 
the generalised coordinates X\ and momenta V\ 



[#RPA j X\] 



-iaj\V\, 



[Hr PAi V\} = iui\X\ 



(8) 



where u>\ are the RPA eigen-frequencies in the rotating 
frame and the associated phonon operators are Ox = 
(Xx - Wx)/V2. Here X\ = E s X sfs, Vx = i£ a P*& 
are bilinear combinations of the quanta ajj~, a& such that 
([fs,g s ']) = V s S SlS > where quantities V s are proportional 
to different combinations of £ f (i = x,+,— ). Further, 
(...) means the averaging over mean field states. Since the 
mean field violates the rotational invariance, among the 
RPA eigen-frequencies there exist two spurious solutions. 
One solution with zero frequency is associated with the 
rotation around the cc-axes, since [H, L x ] = 0. The other 
"spurious" solution at u = Q corresponds to a collec- 
tive rotation, since [ff,£±] = [H : L y ±iL z ] = ^pQL± [p2[. 



The Hamiltonian Eq.([7|) possesses the signature symme- 
try, i.e. [-R^-Hrpa] = (R x = e~ l7ri *), such that it 
decomposes into positive and negative signature terms 



H RPA = H(+)+H(-) 



(9) 



which can be separately diagonalized p2|-|24[ . The nega- 
tive signature Hamiltonian contains the rotational mode 
and the vibrational mode describing the wobbling motion 
p3|,p5[ . We focus on the positive signature Hamiltonian. 
It contains the zero-frequency mode defined by 



[H(- 



—iL x 
Jtv 



,L X ] 



(10) 



and allows one to determine the Thouless-Valatin mo- 
ment of inertia Jtv p6| ]. Here, the angular momentum 
operator L x = J2 S 's/s an d the canonically conjugated 
angle <f> x = i £ s 4>*g s are expressed via f s and g s which 
obey the condition RxdgR" 1 = d s (d s = f s or g s ). Solv- 
ing Eqs.dTq) for the Hamiltonian H(+), 



H(+) = 



k—x,-\-.- 



hu k (ala k + 1/2) -^(Ql + Q[ +)2 



where 



io = 
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(11) 

(12) 
(13) 
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we obtain the expression for the Thouless-Valatin mo- 
ment of inertia 



characterises the collective properties of non-interacting 
fermions, the dynamical moment of inertia reflects the 
changes in the rotating self-consistent mean field due to 
an inter-nucleon interaction. As it was pointed out in 
p7[ , the volume conservation condition, used as a con- 
straint in the mean field calculations, can be interpreted 
as a Hartree approximation applied to an interaction that 
involves the sum of one-, two- etc forces. 



Jtv = Jx 



[2S x qS x 2Sq2 — >5^ (S , 22 



K 2 ' 



(15) 



Here, the term Jx corresponds to the Inglis moment of 
inertia 



Jx = 



E 



E, 



(16) 



The second term in Eq.(|l5|) is a contribution of the 
quadrupole residual interaction in the cranking model. In 
the cranking harmonic oscillator it consists of the terms 
which have the following structure 



Sxm — ^ ^ 



ix my 



n m \r 



0,2 



(17) 



where E s are the energies of particle-hole excitations: 
Ei = 2hcu+, E2 = 2fiuj- 1 E3 = 2hu x , Ei = Hlu+ + Tilu- 
and i?5 = Tilo+ — fouj-. We also introduced the following 
notations: Kn = -t^—k and no = -M-k. 

U 167T A 167T 

The above results are the starting point for our nu- 
merical analysis. To take into account shell effects we 
consider two systems with number of particles A — 
20,64 (N = Z). For = MeV the global mini- 
mum occurs for a prolate shape and for a near oblate 
triaxial shape for A — 20 and 64, respectively pl| . If wc 
trace the configurations which characterise the ground 
states, with increasing rotational frequency both systems 
become oblate. At this point the moment of inertia van- 
ishes, since there is no a kinetic energy associated with 
such a rotation. 

In order to compare various moments of inertia, i.e. 
the Thouless-Valatin, Eq.(|l5|), the Inglis, Eq.([l6|), and 

J$ F = -d 2 E MF /dn 2 with J& A = -d 2 E RPA /dQ 2 
, we calculate the RPA correlation energy E^ A — 
w A — E s ) which includes the positive and neg- 
ative signature contributions. In Figs.[l],|| the results of 
calculations for different moments of inertia are presented 
for A — 20 and 64, respectively. For our knowledge 

this is a first numerical demonstration of the equivalence 

(2) 

between the dynamical moment of inertia J\ 1F calcu- 
lated in the mean field approximation and the Thouless- 
Valatin moment of inertia Jtv calculated in the RPA. 

For the both systems the Inglis moment of inertia Ji is 

(2) 

smaller than the Jtv and J^p and has a different ro- 
tational dependence. While the Inglis moment of inertia 
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FIG. 1. Moments of inertia for TV = Z — 10 system as a 
function of the rotational frequency u> = Q. The definitions 
of different moments of inertia are given in the text. 
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FIG. 2. As in Fig.|l] for N = Z = 32 system 

The sharp drop of all moments of inertia in Fig. 2 is 
caused by the onset of the oblate shape where the col- 
lective rotation does not exist. For A = 64 the onset 
of the oblate deformation occurs at a smaller rotational 
frequency in contrast to the one for the system A = 20. 
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The dynamical moment of inertia Jrpa is larger than 
the Thouless-Valatin moment of inertia. However, from 
our calculations it follows that the contribution of the 
RPA ground state correlations decreases with an increase 
of the number of particles. The difference between the 

(2) 

J R p A and the Jtv is due to the following reason. The 
Inglis moment of inertia is smaller than the Thouless- 

(2) 

Valatin (or J^ F ) value, since the Jtv contains the ef- 
fect of the residual particle-hole interaction. On the other 
hand, the Thouless-Valatin moment of inertia manifests 
the rotational dependence of the residual interaction. 
Thus, we may speculate that inclusion of the phonon 
interaction could help to reproduce the behaviour of the 

(2) 

^ rpa which characterises the rotational dependence of 
the phonon-phonon interaction. 

In summary, using the self-consistent cranking harmonic 
oscillator model, we have numerically proved the equiva- 
lence of the dynamical moment of inertia calculated in 
the mean field approximation to the Thouless-Valatin 
moment of inertia calculated in the RPA. Our result is a 
consequence of the self-consistent condition Eq. which 
minimises the expectation value of the mean field Hamil- 
tonian, Eq. (Q) . This condition is equivalent to the stabil- 
ity condition of collective modes in the RPA ^8|, i.e. uj\ 
to be real and non-negative, and has been used to calcu- 
late different moments of inertia. The rotational depen- 
dence of the both dynamical moments of inertia, J\ lF 

(2) (2) 

and J RP ai is similar, however, the J RFj \ is larger than 

(2) 

the Jmp due to the contribution of the ground state cor- 
relations. This difference between the moments of inertia 
is less important for heavy systems. 
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